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ABSTRACT 


The paper aims to investigate the properties of Brownian motion and Thermophoresis magento hydro dynamics (MHD) stag¬ 
nation point of a Nanofluid Boundary layer flow on a stretching surface with inconsistent thickness. The effect of inconsistent 
thickness is considered and understood that the sheet is defiant. A governing continuity, momentum, angular momentum and 
heat equations mutually with allied boundary conditions are first abridged to a set of self-similar non-linear united ordinary dif¬ 
ferential equations by appropriate transformations. These equations are solved numerically by using the Keller Box method. The 
influence of magnetic parameter M reduced the velocity profile while it increases temperature and nanoparticle volume fraction 
profiles. It is seen that the boundary layer is formed when X>1 and on other hand an inverted boundary layer is formed when <1. 
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INTRODUCTION 

In fluid dynamics, MHD stagnation flow over a stretching 
sheet has many production applications such as tragedy core 
cooling system, glass industrialized, sanitization of crude oil, 
etc. The stagnation-point flow and the flow due to a stretch¬ 
ing sheet are equally essential in theoretical and claim point 
of view. The two-dimensional flow of fluid near a stagnation 
point was first examined by Hiemenz [1]. Later, Chaim [2] 
considering the strain-rate of the stagnation-point flow and 
the stretching rate of the sheet to be identical and found no 
boundary layer structure near the sheet. Whereas, the perfor¬ 
mance of stagnation-point flow over stretching vertical un¬ 
der diverse physical aspects was discussed by Ishak A, Nazar 
R.,[3].Layek et al.,[4], studied the unsteady laminar MHD 
flow and heat Transfer in the stagnation region. The minor¬ 
ity current studies covering similar topics are cited in Refs. 
[5-16]. Besides stagnation point flow, stretching surfaces 


has a ample range of applications in engineering and nu¬ 
merous technical purposes predominantly in metallurgy and 
polymer industry. For instance, plodding cooling of continu¬ 
ous stretching metal or plastic strips can be mentioned which 
have various applications in mass production, pointless to 
say, the final eminence of the product robustly depends on 
the rate of heat transfer from the stretching surface. Crane 
[17] is the first to present a identity-analogous solution in 
the congested analytical form for steady two-dimensional in¬ 
compressible boundary layer flow caused by the stretching 
plate whose velocity varies linearly with the distance from a 
fixed point on the sheet. Following Crane’s study, the ther¬ 
mal approach to this crisis was investigated by Carragher 
and Crane [18]. They tacit that the temperature difference be¬ 
tween the sheets and the ambient is comparative to a power 
of the detachment from the fixed point. Subjecting uniform 
heat flux boundary condition, Dutta et al. [19] presented the 
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temperature allocation for a stretching surface in an ambient 
with different temperatures. 

Nanofluid is a fluid containing small particles, called nano¬ 
particles. These fluids are engineered colloidal suspensions 
of nanoparticles in a base fluidstudied by Buongiomo[20]. 
The nanoparticles used in nanofluids are classically made 
of metals, oxides, carbides, or carbon nanotubes. The base 
fluids used are usually water, ethylene glycol and oilKakac, 
A. Pramuanjaroenkij& Marga, S.J. Palm[21,22]. Recently, 
boundary-layer flow of a nanofluid past a stretching sheet 
was presented by Khan and Pop [23], The natural convective 
boundary-layer flow of a nanofluid past a vertical plate is 
presented systematically by Kuznetsov and Nield [24] Va- 
jravelu et al. [25] presented the exhaustive scrutiny of con¬ 
vective heat transfer in the flow of viscous Ag-water and 
Cu-water nanofluids over a stretching surface. Very recently 
the boundary layers of an uneven stagnation-point flow in 
a nanofluid is investigated by Bachok et al. [26], They an¬ 
alyzed that fluid containing solid particles may drastically 
increase its conductivity. Recently, Bachok et al. coated the 
boundary layers of an unsteady stagnation-point flow in a 
nanofluid. 

Nanofluids have their foremost applications in heat transfer, 
as well as microelectronics, fuel cells, pharmaceutical pro¬ 
cesses, and hybrid-powered engines, domestic refrigerator, 
chiller, nuclear reactor coolant, grinding, space machinery, 
and boiler flue gas temperature reduction. They exhibit en¬ 
hanced thermal conductivity and convective heat transfer 
coefficient counter balanced to the base fluid. Nanotechnol¬ 
ogy has enormous applications in industry. Fluids with nano- 
scaled particles interface are called nanofluid. The term na¬ 
nofluid was proposed by Choi [27] Nanofluid heat transfer 
is a pioneering technology which can be used to augment 
heat transfer. Nanofluid is a deferral of solid nanoparticles 
(l-100nm diameters) inconventional liquids like water, oil, 
and ethylene glycol. The nano particles used in nanofluid are 
usually collected of metals, oxides, carbides, or carbon nano¬ 
tubes. Water, ethylene glycol, and oil are familiar examples 
of base fluids. The effects of the wall thickness constraints 
are velocity power index, ratio of rates of velocities, radiation 
parameter and the Prandtl number on the membrane friction 
coefficient and the local Nusselt number, which represents 
the heat transfer rate at the surface, for the steady stagnation 
point flow and heat transfer towards stretching surface has 
been studied by Prasanna Kumara et al [28] 

In current study we extended the work of Prasanna Kumara 
e.al [28], by taking into consideration the flow of an incom¬ 
pressible nanofluid stagnation point flow towards a stretching 
surface with variable thickness, thermal radiation. The fluid 
determined by a stretching surface located at y = a(x + 
with a fixed stagnation point at x = 0. 



MATHEMATICAL ANALYSES 

Consider the flow of an incompressible viscous 
fluid driven by a stretching surface located at 
y = a(x + b) = with a fixed stagnation point at x = 
0 as shown in figure 1. We imagine that wall is im¬ 
pervious, non-flat with a given profile and the coef¬ 
ficient A being small so that the sheet is adequately 
thin. The stretching velocity UJx) and the ambient 
fluid velocity U(x ) are assumed to be thickness of 
the stretched sheet from the stagnation point i.e., 
UJx ) = u „(x+b) m and U(x) = ujx+b) m , where m is 
the velocity power index. Due to the acceleration 
or deceleration of the sheet, the thickness of the 
stretched sheet may decrease or increase with dis¬ 
tance from the slot, which is reliant on the value of 
the velocity power index. With the above assump¬ 
tions, the boundary layer equations prevailing the 
flow and temperature fields are given by, 
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Where u and v are the velocity components of the fluid along 
x and y directions respectively. M, p,p, k and C p are the mag¬ 
netic parameter, the dynamic viscosity of the fluid, density of 
the fluid, thermal conductivity and specific heat of the fluid 
respectively. T is the temperature, T* is stable temperature of 
the fluid in the in viscid free stream, a is the thermal conduc¬ 
tivity, (pc)p is the effective heat capacity of nanoparticles, 
(pc) f is heat capacity of the base fluid, C is nanoparticle 
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volume fraction, DB is the Brownian diffusion coefficient, 
and DT is the thermophoretic diffusion coefficient. The al¬ 
lied boundary conditions for the present problem are 

“ = y w(4 V = 0. T J X ) = r :o + T q (x + 

m 

c w W= ^ +c a (x + b)r 

L—vTl 

at y = A(x — b)~ ( 5 ) 

u -> u(x ) -,T -> c ->c*, (is y — > go where, = ^o(* + k) m 
is the stretching velocity,UO and b are the physical param¬ 
eter related with stretched surfaces. and denote the 
temperature at the wall and at large distance from the wall 
respectively and TO is the distinctive temperature. is the 
variable wall nanoparticle volume fraction with being a 
constant and is constant nanoparticle volume fraction in 
free stream, To employing the generalized Bernoulli’s equa¬ 
tion, in the free stream U(x) = Ul(x + b)m the equation (2) 
reduces to 


8(d=^rMti = 




where rj is the similarity variable and ^ is the stream funct- 
tion defined as u = and v = - which identically satis¬ 
fies equations (1). Employing the similarity variables (11) 
equations (7) and (10) reduce to the following ordinary dif¬ 
ferential equations 


f"\ri) + mf'Xri) + ~[^~ f’ 2 0 /)] -= 0 ( 12 ) 


+ NbO'WvXtf + NtO’ 2 (rj) = 0 (13) 


0 ”()j) + Le 0 "(»/) = 0 . ( 14 ) 

subject to the boundary conditions (5) which becomes 

f = l 

j 00?) = 1, ft'?) = 1 at i 1 = 0 , 

1/ (v) = A, 8(v) = 0, <p 0?) =0 at JJ (15) 



1 dp 
P d* 
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using (6) into (2) one can obtain 
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Using the Rossel and approximation for radiation, radiation 
heat flux is simplified as 


4c" dr* 
3k* dy 


( 8 ) 


In the above equatio ns, prim e denotes differentiation with 
respect to r\ and P = aJ Jo ( ™ +1: is the wall thickness parameter, 

’ u 

m is the velocity power index, A=zr is the ratio of rates of 

4 _ + t & ye p 

velocities, Nr = " is the radiation parameter and Pr = k 

is the Prandtl number, M= is the magnetic parameter, 
Le = v /°b is the Lewis number. The dimensionless parameter 
Nb (Brownian motion parameter) and Nt (thermophoresis 
parameter) are defined as 


Nb= 


(jMjfc '.c,,. o n,. (pc}p (r„.-r g : 

( ^V v ’ JNt t w {pci f 


where & * and k * are the Stefan-Boltzman constant and mean 
absorption co-efficient, respectively. Assuming that the tem¬ 
perature differences within the flow such that the term T 4 
may be expressed as a linear function of the temperature, 
we expand in a Taylor series about T.^ and neglecting the 
higher order terms beyond the first degree in (T-T :ci ) we get 

T 4 G= 4TjT - 3Tj* ( 9 ) 


Substituting equations (8) and (9) in equation (3) reduces to 


pc p 



IS tr*T^ d 2 T 
3 k* dy” 


( 10 ) 


Here Pr, Le, Nb, Nt, M,and Nr denote the Prandtl number, the 
Lewis number, the Brownian motion parameter, the thermo¬ 
phoresis parameter, magnetic parameter, Radiation param¬ 
eter respectively. This boundary value problem is reduced to 
the classical problem of flow and heat and mass transfer due 
to a stretching surface in a viscous fluid when n = 1 and Nb 
= Nt = 0 in eqs (13) and (14). 

The quantities of practical interest, in this study, are the local 
skin friction C, Nusselt number N and the Sherwood num- 

fx ux 

ber Sh which are defined as 




pul 


7 — ,Nu 


y =o 


xq w 


k(T w ~ T J 


, Sh 


xq m 


D b (C w -CJ 


(16) 


The momentum and energy equations can be transformed us¬ 
ing the following correspondence renovation 
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Where k is the thermal conductivity of the nanofluid, and 
q w , Vm are the heat and mass fluxes at the surface, respec¬ 
tively, given by 


q w = - 
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Substituting Eq (11) into Eqs (16) - (17), we obtain 
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R e~/ 2 Nu x 



Re -Sh = - 


n +1 


^( 0 ) 


Where Re x = u w ~ is the local Reynolds number. 


Table 1: Comparison of the values of membrane fric¬ 
tion coefficient / "( 0 ) for a choice of values of m and 
for fixed values X= Pr= Nr = o and = OS 


M 

Fang et al. 

Prasanna 
Kumar et al 

Present results 

10 

-1.0603 

-1.06034 

-1.0602 

7 

-1.0550 

-1.05506 

-1.0549 

3 

-1.0359 

-1.03588 

-1.0358 

1 

-1.0000 

-1.00000 

-1.0000 


Table 2: Result Values of skin friction -/"(0), Nusslet 
number -G ( 0 ) and Sherwood number - 0 '( 0) for differ¬ 
ent values of physical parameters 


£ 

Pr 

Nr 

M 

Nt 

Nb 

Le 

A 

- no ) 

-£"(o) 

-0'(O) 

0.5 

3 

2 

2 

1 

1 

2 

O.5 

0.^022 

O.4178 

I.414O 

1.0 








0.6620 

0*3933 

1.2632 

i -5 








0.6241 

0.3707 

I.1254 

0.5 

1 







0.7022 

0.3481 

I.414O 


2 







0.7022 

0.4043 

1.4298 


3 







0.7022 

0.4178 

I.4449 



2 






0.7022 

0.4178 

I.414O 



2-5 






0.7022 

0.4144 

I.4072 



3 






0.7022 

0.4091 

1.4028 




1 





0.6676 

0.3990 

I.6180 




2 





0.7022 

0.4178 

I.414O 




4 





O.7272 

0.4340 

1.2630 





1 




0.7022 

0.4178 

I.414O 





2 




0.7022 

0.3501 

1.4308 





3 




0.7022 

0.3040 

1.4562 






O.5 



0.7022 

0.5146 

I.261O 






1 



0.7022 

0.4178 

I.414O 






i -5 



0.7022 

0.3405 

1.4546 







1 


0.7022 

0.4574 

O.9255 







2 


0.7022 

0.4178 

I.414O 







3 


0.7022 

0.3979 

I.7520 








0.2 

0.9528 

0.3874 

I.3OIO 








1-5 

O.9797 

0.5176 

1.7763 








2 

2.1870 

0.5629 

I.94OI 


RESULTS AND DISCUSSION 

The classification of equations (l)-(4) along with the bound¬ 
ary condition (5) has been solved numerically by using im¬ 
plied finite difference scheme known as Keller Box method. 

The abridged Eqs. (12)—(15) are nonlinear and united, and 
thus their exact methodical solutions are not potential. They 
can be solved numerically using Keller - Box for dissimilar 
values of parameters such as magnetic parameter, Prandtl, 
Eckert and the Lewis numbers, the Brownian motion pa¬ 
rameter and the thermophoresis parameter. The possessions 
of the emerging parameters on the dimensionless velocity, 
temperature, membrane friction, the rates of heat and mass 
transfer are investigated. 

The important steps in using the Keller Box method are: 

1) Tumbling higher order ODEs (systems of ODEs) in to 
system of first order ODEs; 

2) Writing the systems of first order ODEs into differ¬ 
ence equations using central differencing scheme; 

3) Liberalizing the difference equations using Newton’s 
method and writing it in vector form; 

4) Solving the system of equations using block elimina¬ 
tions method. 

In order to solve the above degree of difference equations 
numerically, we adopt Matlab software which is very effi¬ 
cient in using the well known Keller Box method. 

In order to study precision of our results compared the pre¬ 
sent results with those of PrasannaKumaraet.al.,[28],when 
neglect the effects of Nt and Nb. The evaluation results en¬ 
sures there is a good concurrence between the present and 
previous works Prasanna Kumara et al [28]. The computed 
results for some values of the prevailing flow parameters 
M,m, p y A Nr, Nt, Nb ,Le are analyzed graphically for veloc¬ 
ity, temperature and nanoparticle volume fraction profiles in 
figures (2) to (9). 

Fig. 2(a-c) exhibits the effect of magnetic parameter on the 
velocity. It is empirical that the velocity profile of the fluid 
is inconsequentially reduced with increasing values of M. 
An increase in magnetic parameter M results in a strong 
diminution in velocity. This is due to the fact that magnetic 
field introduces a retarding body force which acts oblique 
to the direction of the applied magnetic field. This body 
force, known as the Lorentz force, decelerates the boundary 
layer flow and the boundary layer thickness decreases with 
increase in magnetic parameter. The authority of magnetic 
parameter M on the temperature and concentration profiles is 
clearly observed to be appreciably enhanced with increasing 
magnetic parameter. As the Lorentz force is a resistive force 
which opposes the fluid motion, so heat is produced and as a 
result, the thermal boundary layer thickness become thicker 
for stronger magnetic field. 
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Fig 3(a)-(c), represents the deviation of velocity, temperature 
and nanoparticle volume fraction for different values of m. It 
is found from fig 3(a) that the velocity profiles decrease with 
the increase of m in case of A < 1 and increase the velocity 
profiles with increase of m in case of ^ > L For a fixed value 
of A = 0.5 the temperature and nanoparticle volume fraction 
increases with the increase of m is depicted from the figure. 
One can seen that if m= 1, the problem reduces to a flat sheet 
problem. 

Fig 4(a)-(c), illustrates the deviation in the velocity, tempera¬ 
ture and nanoparticle volume fraction, respectively, for the 
different values of P. From the figure it is apparent that if P 
increases the fluid velocity is increases for a fixed value of 
A < 1 whereas the fluid velocity decrease with increase of P 
for the value of (A > i). This is because for higher value of P 
the boundary layer becomes thicker. It can be experimental 
from the figure that the temperature profiles and nanoparticle 
volume fraction profiles increases with increase of P. It can 
be noticed that an increase in wall thickness parameter re¬ 
sults, increases the profiles of temperature and nanoparticle 
volume fraction. 

Fig 5(a)-(c), exhibits the ratio of frees tiring velocity param¬ 
eter A on velocity, temperature and nanoparticle volume frac¬ 
tion profile respectively. It is observed that there is a decrease 
in the velocity and nanoparticle volume fraction profiles is 
noticed with an increasing A. It is originated that when the 
stretching velocity is less than the free string velocity i.e., 
A < 1, the flow as a boundary layer structure actually means 
the straining motion near the stagnation region increases, so 
the acceleration of the exterior stream increases which leads 
to decrease in the width of the boundary layer with increase 
in A. In case of A < 1 i.e., the stretching velocity of the float 
up exceeds the free string velocity then overturned boundary 
layer structure is formed. For ^ 1 there is no boundary 

layer formation because the stretching velocity is equal to 
free string velocity. It is apparent from the figure that the 
boundary layer thickness decreases with an increase in A. 

Fig 6(a)-(b), predicts the influence of radiation parameter Nr 
on temperature and nanoparticle volume fraction profiles, re¬ 
spectively. The radiation parameter Nr defines the compara¬ 
tive contribution of conduction, heat transfer to thermal ra¬ 
diation transfer. It is evident that an increase in the radiation 
parameter results with an increase in temperature field. It is 
seen that the nano particle volume fraction profiles decreases 
with the increase of radiation parameter Nr. This is an agree¬ 
ment with physical fact the thermal boundary layer thickness 
increases with an increase in Nr. 


Fig7(a)-(b), illustrates the dimensionless temperature and 
nanoparticle volume fraction profiles for diverse values of 
thermophoresis parameter Nt. It is found from fig 7(a) the 
temperature increases with increase in nanoparticle volume 
fraction Nt. It is observed from fig 7(b) that increasing the 
value of thermophoresis Nt leads to enhances the nanoparti¬ 
cle volume fraction profiles. 

Fig 8(a)-(b), has been plotted to displayed the Brownian 
motion parameter Nb on temperature, nanoparticle volume 
fraction profiles, respectively. It illuminates that the temper¬ 
ature profile increases with an increasing Brownian motion 
parameter Nb. It is also noticed from 8(b) that the nanopar¬ 
ticle volume fraction profiles decrease of increase Brown¬ 
ian motion parameter Nb. The concentration boundary layer 
reduces with an increasing Nb. 

Fig 9(a)-(b), depicts the effects of Lewis number Le on tem¬ 
perature field and nanoparticle volume fraction profiles re¬ 
spectively. It is seen that increasing values of Lewis number 
Le result decrease the temperature field along with thermal 
boundary layer thickness. It is also observed that the high 
end profiles are to decrease with increasing Le. Thus increas¬ 
ing Le results to decrease the diffusion, which finally yields 
into a decrease of nanoparticle volume fraction. 



Figure 2: Effects of M on Velocity, Temperature and Nanopar¬ 
ticle volume. 
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Figure 3: Effects of m on Velocity, Temperature and Nanopar¬ 
ticle volume. 



T1 (a) 




T! (c) 


Figure 4: Effects of/? on Velocity, Temperature and Nanopar¬ 
ticle volume. 



(a) 




Figure 5: Effects of A on Velocity, Temperature and Nanopar¬ 
ticle volume. 




T1 (b) 

Figure 6 : Effects of Nr on Temperature and Nanoparticle vol¬ 
ume. 
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T1 ( a ) 



n (b) 

Figure 7: Effects of Nt on Temperature and Nanoparticle volume. 




■n ( b ) 

Figure 8: Effects of Nb on temperature and Nanoparticle volume. 




Figure 9: Effects of Le on Temperature and Nanoparticle volume. 


CONCLUSION 

The sturdy of two-dimensional boundary layer flow of an 
incompressible viscous nanofluid determined by a stretching 
surface with a fixed MHD stagnation point are present. The 
consequential partial differential equations are condensed to 
ordinary differential equations by using appropriate compar¬ 
ison transformations. The classification of ordinary differen¬ 
tial equations are solved numerically by Keller Box method 
the numerical results obtainable to analyze the flow and heat 
transfer distinctiveness, Some of the imperative finding of 
four investigations are, 

1. The effects of M, m and leads to enhance the pro¬ 
files of Temperature and Nanoparticle volume frac¬ 
tion, while the Magnetic parameter decreases velocity 
profile. 

2. The stretching Velocity is less than the free stream 
velocity A > 1. The flow as a boundary layer structure 
actually tells the straining motion near the stagnation 
region. 

3. The stretching Velocity of the surface exceeds the 
free stream velocity A < l transformed boundary layer 
structure is formed. For A =1 there is no boundary lay¬ 
er configuration. With the increase of 2, Temperature 
and Nanoparticle volume profiles decreases. 

4. The effects of Radiation parameter (Nr) enhance the 
Temperature profiles and Nanoparticle volume frac¬ 
tion profiles. 

5. The Temperature profiles and Nanoparticle volume 
fraction profiles increases with the effect of Thermo¬ 
phoresis parameter. 

6. The augmentation of Brownian motion parameter is to 
develop the temperature profiles and reduce the Nano¬ 
particle volume fraction. 

7. The f rr (0) value decreases with increase of P and m. 
But there is no changes when increases the values of 
Nr, Nt, Pr changes. The -0 f (o) value decrease with in¬ 
crease in the values Nb, Nt, and Le. 
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